Astronomy & Astrophysics manuscript no. 7015'galvez'fV 


© ESO 2008 


February 1, 2008 





o 
o 

(N 



Multiwavelength optical observations of chromospherically active 

binary systems 

V. FF UMa (2RE J0933+624): a system with orbital period variation *,' 

M.C. Galvez 1 ' 2 , D. Montes 1 , M.J. Fernandez-Fi gueroa 1 , E. De Castro 1 , and M. Cornide 1 



1 Departamento de Astrofi'sica, Facultad de Ciencias Ffsicas, Universidad Complutense de Madrid, E-28040 Madrid, 
Spain 

E-mail: mcz@astrax.fis.ucm.es 

2 Department of Astronomy, University of Florida, Gainesville, FL 32611, USA 



3 

in 



Received .... 2006 / Accepted .... 2007 

ABSTRACT 



Context. This is the fifth paper in a series aimed at studying the chromospheres of active binary systems using several 
optical spectroscopic indicators to obtain or improve orbital solution and fundamental stellar parameters. 
Aims. We present here the study of FF UMa (2RE J0933+624), a recently discovered, X-ray/EUV selected, active 
binary with strong Ha emission. The objectives of this work are, to find orbital solutions and define stellar parameters 
from precise radial velocities and carry out an extensive study of the optical indicators of chromospheric activity. 
Methods. We obtained high resolution echelle spectroscopic observations during five observing runs from 1998 to 2004. 
We found radial velocities by cross correlation with radial velocity standard stars to achieve the best orbital solution. 
We also measured rotational velocity by cross-correlation techniques and have studied the kinematic by galactic space- 
velocity components (U, V, W) and Eggen criteria. Finally, we have determined the chromospheric contribution in 
optical spectroscopic indicators, from Ca n H & K to Ca II IRT lines, using the spectral subtraction technique. 
Results. We have found that this system presents an orbital period variation, higher than previously detected in other 
ly-^ ■ RS CVn systems. We determined an improved orbital solution, finding a circular orbit with a period of 3.274 days. 

\^ ' We derived the stellar parameters, confirming the subgiant nature of the primary component (JVfp = 1.67 Mq and 

\^») , Rsinip = 2.17 Rq) and obtained rotational velocities (vsini), of 33.57±0.45 km s" 1 and 32.38±0.75 km s" 1 for the 

. primary and secondary components respectively. From our kinematic study, we can deduce its membership to the Castor 

moving group. Finally, the activity study has given us a better understanding of the possible mechanisms that produce 
the orbital period variation. 

o ■ 

' Key words, stars: FF UMa - stars: 2RE J0933+624 - stars: activity - stars: binaries:: spectroscopic - stars: chromo- 

f"*^ , spheres - stars: late-type 

> 

■ 1. Introduction provide the information for several optical spectroscopic 

5h ' . . . . ... features that are formed at different heights in the chro- 

d • This paper is a continuation of our ongoing project aimed at here (see Montes et aL 1997 p I; Montes et al. 

studying the chromospheres of active binary systems using igg8 Paper n Monteg ^ ^ 2m p&per m Galvez et 

multiwavelength optica l observations. These observations aL 2QQ ^ p aper IV) In addition to studying ste n ar activ _ 

Send offprint requests to: M.C. Galvez our hi g h resolution spectroscopic observations allow us 

* Based on observations collected with the 2.2 m telescope at to determine radial velocities and obtain and improve fun- 

the Centro Astronomico Hispano Aleman (CAHA) at Calar Alto damental stellar parameters. While several systems have 

(Almen'a, Spain), operated jointly by the Max-Planck Institut been studied, this is the first time we have found an orbital 

fur Astronomie and the Instituto de Astrofi'sica de Andalucfa period variation, giving us new clues into the understand- 

(CSIC); with the Nordic Optical Telescope (NOT), operated on i ng f activity-orbit relation. When combined with other 

the island of La Palma jointly by Denmark, Finland, Iceland, examples, the study of this type of system could help us 

Norway and Sweden, in the Spanish Observatorio del Roque de un d er stand how orbital dynamics are affected by physical 

Los Muchachos of the Instituto de Astroffaica de Canarias; with sses intrinsic to the w tem (Lanza 2006) 

the 2.1 m Otto Struve Telescope at McDonald Observatory of v ; 
the University of Texas at Austin (USA) and with Hobby-Eberly 
Telescope, which is a joint project of the University of Texas at 

Austin, the Pennsylvania State University, Stanford University, We focus here on the X-ray/EUV selected chromospher- 

Ludwig-Maximilians-Universitat Miinchen, and Georg-August- ically active binary FF UMa (2RE J0933+624, HD 82286, 

Universitat Gottingen. SAO 14919). It is an SB2 system with V = 8.35 mag. 

Tables 8 and 9 are only available in electronic form via First classified by Jeffries et al. (1995) as two G5V or 



http://www.edpsciences.org G5V/G5IV stars, it was reclassified by Henry et al. (1995) 



Correspondence to: mcz@astrax.fis.ucm.es and Strassmeier et al. (2000) as an K0IV/K0IV. 
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Table 1. Observing log 
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Henry et al. (1995) reported a photometric period of 
3.270 days obtained from a periodogram analysis derived 
from photometric data, and estimated a rotational velocity 
vs'mi = 35 km s _1 for both components. Jeffries et al. 
(1995) obtained an orbital period of 3.28 days from 15 radial 
velocity measurements and suggested an eccentricity less 
than 0.18. Their estimated value of i?sini was in agreement 
with a subgiant primary. Strassmeier et al. (2000) found a 
rotational period of 3.207 days and photometric variation 
amplitudes of AU « 0.15 mag. All previous authors have 
reported that this star presents very strong chromosphcric 
activity and the Ha emission line is detected above the 
continuum for both components. 

In this paper, we present high-resolution echelle spec- 
tra of this system. We measured radial velocities using the 
cross-correlation technique and obtained an orbital period 
variation during 11 years of observations. In spite of this 
variation, we achieved a good orbital solution, finding that 
the mean orbital period is similar to the photometric one, 
indicating synchronous rotation. 

In addition, we applied the spectral subtraction tech- 
nique to study the chromospheric excess emission in the 
Ca II H & K, Ca II IRT, Ha and other Balmcr lines 
of the primary and secondary components of the system. 
Preliminary results for this system can be found in Galvez 
(2005); and Galvez et al. (2006, 2007). 

In Sect. 2 we give the details of our observations and 
data reduction. In Sect. 3 we discuss the nature of the or- 
bital period variation and give the orbital solution of the 
binary system. In Sect. 4 we give the derived stellar and 
kinematic parameters. The behavior of the different chro- 
mospheric activity indicators is described in Sect. 5. Finally, 
in Sect. 6 we present our conclusions. 

2. Observations and data reduction 

We obtained high resolution echelle spectra of FF UMa 
during five observing runs from 1998 to 2004: 

1) 2.1 m-SANDIFORD, McDonald Obs., 1998/01 
During this observing run, which extended from 12 to 21 
January 1998, we used the 2.1 m Otto Struve Telescope 
at McDonald Observatory Texas (USA) with the Sandiford 
Cassegrain Echelle Spectrometer (SCES), equipped with a 
1200x400 pixel CCD detector. The wavelength range cov- 
ers from 6400 to 8800 A in 31 orders. The reciprocal dis- 



persion ranges from 0.06 to 0.08 A/pixel and the spectral 
resolution, determined as the full width at half maximum 
(FWHM) of the arc comparison lines, from 0.13 to 0.20 A. 
In the fourth night, the central wavelength was changed to 
include the Na I D b D 2 (5889.95, 5895.92 A) and He I D 3 
(5876 A) lines. Therefore the wavelength range changed to 
5600-7000 A. 

2) 9.2 m-HET, McDonald Obs., 2000/01 

We used the 9.2 m Hobby- Eberly Telescope (HET) and the 
medium resolution spectrograph UFOE (Upgraded Fiber 
Optic Echelle) equipped with a 1200x400 pixel CCD de- 
tector, located at McDonald Observatory Texas (USA) on 
22 to 24 January 2000. The wavelength range covers from 
4400 to 9150 A in 26 orders. The reciprocal dispersion 
ranges from 0.06 to 0.17 A/pixel and the spectral resolution 
(FWHM) ranges from 0.14 to 0.42 A. 

3) and 4) 2.2 m-FOCES, CAHA, 2002/04 and 2004/04 
We utilized the Fibre Optics Cassegrain Echelle 
Spectrograph (FOCES) (Pfciffer et al. 1998) with a 
2048x2048 24^, SITE#ld CCD detector on the 2.2 m tele- 
scope at the German Spanish Astronomical Observatory 
(CAHA) (Almeria, Spain) to obtain spectra between 22 
to 26 April 2002 and from 29 March to 7 April 2004. The 
wavelength range covers from 3450 to 10700 A in 112 
orders. The reciprocal dispersion ranges from 0.04 to 0.13 
A/pixel and the spectral resolution (FWHM) ranges from 
0.08 to 0.35 A. 

5) 2.56 m-NOT-SOFIN, Roque dc los Muchachos Obs., 
2004/04 

We used the 2.56 m Nordic Optical Telescope (NOT) lo- 
cated at the Observatorio del Roque de los Muchachos (La 
Palma, Spain) on 2 to 6 April 2004. We used The Soviet 
Finnish High Resolution Echelle Spectrograph (SOFIN) 
with an echelle grating (79 grooves/mm), ASTRO MED- 
3200 camera and a 2052x2052 pixel 2K3EB PISKUNOV1 
CCD detector. The wavelength range covered from 3545 to 
10120 A in 42 orders. The reciprocal dispersion ranges from 
0.033 to 0.11 A/pixel and the spectral resolution (FWHM) 
from 0.14 to 0.32 A. We note that we had some problems 
with the wavelength calibration, during the arc lamp spec- 
tra exposures, and as a consequence, could not rely on the 
absolute wavelength calibration of the spectra taken in this 
observing run. Therefore, these spectra have not been used 
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Table 2. Radial velocities 
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115 


-4.99 ± 6.20 


- 


HET00 


51561.943 
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FOCES04 


53102.5640 


93 


11.02 ± 3.46 


-40.39 ± 8.89 



1 JEF(95): Jeffries et al. (1995) 

to determine radial velocities, although we have used them 
for the remaining analysis. 

In Table Q] we present the observing log. For each obser- 
vation we list date, UT, and the signal to noise ratio (S/N) 
obtained in the Ha line region. 

We extracted spectra using the standard reduction pro- 
cedures in the IRAfQ package (bias subtraction, flat-field 
division and optimal extraction of the spectra) . The wave- 
length calibration was obtained by taking spectra of a Th- 
Ar lamp. Finally, we normalized the spectra by a polyno- 
mial fit to the observed continuum. 



3. Orbital period variation 

3.1. Radial velocities 

We determined the heliocentric radial velocities by making 
use of cross-correlation technique (see Paper IV) . The spec- 



1 IRAF is distributed by the National Optical Observatory, 
which is operated by the Association of Universities for Research 
in Astronomy, Inc., under contract with the National Science 
Foundation. 
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Fig. 1. CCF of FF UMa (2RE J0933+624) in FOCES04 
observing run. Irregular profiles can be seen in the two 
peaks. These irregularities can produce significative errors 
in radial velocity determination. 



tra of the target were cross-correlated order by order, using 
the routine FXCOR in IRAF, against spectra of radial veloc- 
ity standards with similar spectral type taken from Beavers 
et al. (1979). We derived the radial velocity for each order 
from the position of peak of the cross-correlation function 
(CCF) and calculated the uncertainties based on the fit- 
ted peak height and the antisymmetric noise as described 
by Tonry & Davis (1979). As FF UMa is an SB2 system 
we note two peaks in the CCF, associated with the two 
components, and fit each one separately. When the com- 
ponents are too close, we used deblending fits. It is worth 
mentioning that the uncertainties returned by FXCOR for 
SB2 binaries are overestimated; when fitting each star, the 
presence of the other will increase the antisymmetric noise, 
thereby biasing the error. 

As Fig. Q] shows, the irregular profiles of the CCF (dou- 
ble peaks and asymmetries) can produce significant errors 
in radial velocity measures. These irregularities may come 
from photospheric activity features on the stellar surface 
of both components that disturb the profile of the pho- 
tospheric lines and induce variations in the peak of the 
CCF. However, this behavior may also be due to the dif- 
ference in rotational velocity (wsini) between the problem 
and radial velocity star. When the spectrum of the stan- 
dard star was broadened to the same rotational velocity 
of FF UMa (vsini « 30 km s" 1 ) the profiles of CCF be- 
come smoother and could be fit with a Gaussian profile, 
see Fig. [5J Therefore all the radial velocities given in this 
paper have been calculated by cross-correlation with this 
rotational broadened spectrum of the standard star. 

In Table [5] we list, for each spectrum, the heliocentric 
radial velocities (Vhel) and their associated errors (try) ob- 
tained as weighted means of the individual values deduced 
for each order in the spectra. Those orders which contain 
chromospheric features and prominent telluric lines have 
been excluded when determining the mean velocity. 
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Fig. 2. CCF of FF UMa (2RE J0933+624) in FOCES04 
observing run obtained when we broad the standard star 
to FF UMa rotational velocity. Irregular profiles become 
smoother and can be fitted with a Gaussian. 



Table 3. Variations of Orbital Period (P) 



Year 


Tconj 


O-C 


dP/P 




HJD 


day 






(2400000 +) 






1993 


49055.8789 


-2.0089 


4.979x10-4 


1998 


50824.4023 


-1.1755 


5.187x10-4 


2000 


51561.2227 


-0.9005 


5.887x10-4 


2002 


52386.6758 


-0.3615 


5.134x10-4 


2004 


53090.8398 







3.2. T con j variations 

With 35 radial velocity data from our measures and nine 
from Jeffries et al. (1995) (see Table [5]), we computed the 
orbital solution of this system. Although we obtained good 
results when we fit orbital solution for each observing run 
data separately, some orbital parameters changed from one 
epoch to another. When we tried to fit the orbital solution 
with all the data, we could not find any satisfactory result. 

We decided to recalculate the orbital solution of each ob- 
serving run, using the period obtained from the FOCES04 
observing run and assuming a circular orbit (since e is only 
~ 10~ 2 ). We determined that the solutions are very similar 
except at time of conjunction, T con j, (see Fig.|3J. Therefore 
we shifted in phase the solutions of every run taking as stan- 
dard the FOCES04 solution. As we can see in Fig.H all data 
points are now in agreement with the orbital solution. We 
obtained the phase shift calculating the conjunction time 
differences between the conjunction time obtained in ev- 
ery run fit and the conjunction time of FOCES04 run fit, 
(O — C) (Observed - Calculated = T con j difference coming 
from orbital solution fit in each run and the FOCES04 run). 

In Fig. O we represent the temporal variation of the 
O — C (T con j) for every run. The maximum amplitude of the 
variation found in our data is 2 days and shows a decreasing 
linear tendency. We would need a longer temporal range 



of observations to test if the tendency remains linear or 
becomes sinusoidal as one might expect if there is a cycled 
behavior related with the activity cycle (see Frasca & Lanza 
2005 and Sect. 3.3). 

When the O — C (T con j) variations are transformed to 
relative orbital period variations, we find dP/P m 10~ 4 in 
11 years, that is, one order of magnitude higher than the 
largest value observed until now in HR 1099, (see Table[3|). 

3.3. Discussion 

To explain the behavior described above, we considered sev- 
eral options: 

— First, we explore how the existence of a third distant 
star as a component of the system could modify the 
main orbit. In Fig. we plotted the center of mass 
radial velocity, 7, obtained for each observing run, ver- 
sus time. The amplitude of variations in 7 amounted 
to 3 km s -1 , over 11 years. Such differences are large 
enough that they are unlikely to be due to zero-point 
(instrumental) differences between different runs. The 
variations in 7 could indicate the presence of a third 
component; if the third star is small and its period long, 
the reflex motion of the binary will necessarily be small. 
As an example, using eq.(30) from Cumming (2004), a 
third star with a mass between 0.4 and 0.6 Mq and an 
orbital period of 20 to 40 yr in an edge-on circular orbit 
would produce a semi-amplitude K on the binary be- 
tween 2.1 and 3.8 km s" 1 . In addition, if the orbit were 
significantly eccentric, as is often the case for such long 
periods, the K amplitude could be larger. The presence 
of the third component could easily induce a change in 
7 similar to that observed. Therefore with the present 
data we cannot dismiss the possibility that these varia- 
tions are due to a third body. 

— Another explanation of our observations could be an 
orbital modulation due to the variation of activity with 
time, explained as a consequence of cyclical variations 
of the quadrupole-moment of both components of the 
system during the magnetic activity cycle. This mech- 
anism presented by Matese & Whirtmere (1983) and 
developed by Applegate (1992) and Lanza et al. (1998), 
has been used in the study of several RS CVn systems 
such as SZ Psc (Kalimeris et al. 1995), RT Lac, RS CVn, 
WW Dra, etc. (Lanza & Rodono 1999), XY UMa 
(Sowell et al. 2001) and HR 1099 (Garcfa-Alvarez et 
al. 2003; Frasca & Lanza 2005, Lanza 2006). 
Applegate (1992), described the initial model in which 
the orbital period variation is due to the gravitational 
coupling of the orbit to changes in the quadrupole mo- 
ment (rotational oblateness) of a magnetically active 
star in the system. The quadrupole moment of a star 
is determined by the rotation rate of its outer layers -if 
angular momentum is transferred to the outer layers, 
they rotate faster and the star becomes more oblate. 
The gravitational acceleration varies if the shape varies; 
this shape variation is measured by the change of the 
quadrupole moment of the star. On the contrary, if 
the outer layers loses angular momentum, the oblate- 
ness decreases. As the dynamo mechanism implies the 
qualitative shearing of magnetic field by differential ro- 
tation, the last should vary through the activity cy- 
cle. Applegate (1992) says that quantitatively, a sub- 
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Orbit of RE J0933 (FF UMa) 
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Fig. 3. Example of the phase shift between orbital solution of FOCES04 and McDonald98 observing runs. Radial velocities 
and orbital fit (solid line) of FOCES04 is plotted and radial velocities of McDonald98 are superimposed. Filled symbols 
correspond to primary and open symbols to secondary. 



Orbit of RE J0933 (FF UMa) 
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Orbital Phase (T conj ) 

Fig. 4. Radial velocities of all observing runs. The orbital solution fit of FOCES04 observing run is plotted here (solid 
line) and the radial velocities of the rest runs (shifted in phase) are superimposed. See text for explanation (Sect. 3.2). 



surface magnetic field of several kilogauss can exert a 
large enough torque to transfer the angular momentum 
needed to make the observed period changes. 
Lanza et al. (1998) studied several possibilities of the 
process to explain the period variations in RS CVn 
systems with different kind of dynamos. They showed 
that variations of about 100 Gauss in a poloidal mag- 
netic field could produce the observed variations, while 
Applegate suggested variations in order of kilogauss. 
Lanza & Rodono (1999) compiled 46 binary systems of 
different types (RS CVn, WW UMa, etc.) to evaluate 
the effects of the quadrupole moments change. 
Lanza (2005), analyzes the Applegate model predictions 
and the observed results in RS CVn stars. He suggested 
that the model should be rejected because it fails to ex- 
plain the orbital period variations of classical RS CVn 
close binary systems. The required variation of the inter- 
nal differential rotation is too large to both agree with 
the observations and oppose turbulent dissipation. He 



concludes that any similar hypothesis to explain this 
phenomenon should include the effect of the Lorentz 
force on the gravitational quadrupole moment, or, that 
an entirely new theoretical framework is needed to in- 
terpret the observed orbital period variations in mag- 
netically RS CVn binaries. 

Based on the Lanza (2005) review of the Applegate 
model, Frasca & Lanza (2005) and Lanza (2006), con- 
tinued with the characterization of the orbital period 
variation of HR 1099. They suggested that there is an 
interaction between the magnetic fields of the Kl IV 
subgiant (primary component), and the magnetic fields 
of the G5 V component (secondary). In the primary, 
the hydromagnetic dynamo action is maintained in the 
deep fast-rotating convective envelope, while in the sec- 
ondary, the magnetic field comes from an outer con- 
vection zone with a smaller radial extension, implying 
that its dynamo is less efficient (reflected by its lower 
level of activity). They argued that the Applegate clas- 
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sic model could not explain the variation found while 
their assumptions of the relation between orbital pe- 
riod cycle and the activity cycle could explain the large 
variation measures. They mention, however, the need 
to verify their claim with a larger study including other 
systems. 

Summing up, the results found until now indicate that 
variations in the orbital period based on the Applegate 
model should be revised and that the strong temporal 
activity observable changes could reflect the relation be- 
tween the orbital period variations and the changes in 
magnetic field and gravitational quadrupolar moment. 
The variations found in previous RS CVn systems are 
about dP/P w 10~ 6 — 10~ 5 in a 7 to 109 years range. 
The high level of chromospheric activity of both com- 
ponents of FF UMa (spectral types Kl IV and K0 V, 
see Sects. 4.2 and 5) could imply a strong interaction 
between larger and more efficient dynamos. This could 
explain the order of magnitude difference between the 
orbital period variation of this system and the one de- 
tected in other RS CVn systems like HR 1099. 
— Finally, although the above explained activity-related 
period variation is our preferred interpretation, we 
should mention that in the case of the presence of a 
third component, there could well be changes in the el- 
ements of the inner orbit. If the eccentricity were small 
but non-zero there could be apsidal motion in the binary 
due to the third body and the change in the longitude 
of periastron would be seen as a change in the time of 
conjunction. 

3,4. Orbital solution 

As a consequence of the results in Sects. 3.2 and 3.3, we have 
computed the orbital solution of this system using radial ve- 
locity data from the FOCES04 observing run. We chose this 
run because it has a large number of data points (13) and 
superior spectral resolution. The radial velocity data are 
plotted in Fig. |4j Solid symbols represent the primary and 
open symbols represent the secondary. Each observing run 
is represented with a different symbol. The curve represent 
a minimum \ 2 fit orbit solution. The orbit fitting code uses 
the Numerical Recipes (Press et al. 1986) implementation 
of the Levenberg-Marquardt method of fitting a non-linear 
function to the data, which weights each datum according 
to its associated uncertainty. The program simultaneously 
solves for the orbital period, P rb, the epoch of periastron 
passage, T con j, the longitude of periastron, u, the eccentric- 
ity, e, the primary star's radial velocity amplitude, Kp, the 
heliocentric center of mass velocity, 7, and the mass ratio, 
q. The orbital solution and relevant derived quantities are 
given in Table|U In this table, we give T con j as the heliocen- 
tric Julian date of conjunction with the hotter star behind 
the cooler star, in order to adopt the same criteria used 
in previous papers. We used this criterion to calculate the 
orbital phases of all the observations reported in this paper. 

This binary results in a circular orbit (adopted) with 
an orbital period of about 3.274 days. Since P p hot ~ 3.270 
days, we can say that it is a synchronous system. The mass 
ratio of 2.12 calculated let us conclude that the compo- 
nents have a different spectral type. The obtained param- 
eters are in agreement with the values reported by Jeffries 
et al. (1995). 



Table 4. Orbital solution 



Element 


Value 


Uncertainty 


Units 


Porb 


3.274 


0.054 


days 


-^conj 


53090.84 


0.18 


HJD (2400000 +) 


LU 


0.00 


0.00 


degrees 


e 


0.00 


0.00 


(adopted) 


K P 


29.55 


0.95 


km s 


K s 


62.52 


3.60 


km s 


1 


-3.23 


0.77 


km s — 


q = Mp/Ms 


2.12 


0.10 




dp sim 


1.330 


0.048 


10 6 km 


IMj iJ±±LL 


2.81 


0.17 


10 6 km 


a sim 


4.14 


0.18 


10 6 km 


if 


0.028 




AU 


r> 


5.95 




Rq 


M P sin 3 i 


0.180 


0.024 


M 


Ms sin 3 i 


0.085 


0.012 




f(M) 


0.00875 


0.00099 


Mq 



4. Stellar Parameters of the binary system 

We give the adopted stellar parameters of FF UMa in 
Tabled The photometric data (B - V, V, P pho t) are taken 
from SIMBAD, Jeffries et al. (1995), Henry et al. (1995) 
and Strassmeier et al. (2000). Orbital period (P rb) aud 
projected rotational velocity (v sin i) have been determined 
in this paper (see Sects. 3.4 and 4.1). The astrometric data 
(parallax, 7r; proper motions, [i a cos6 and fig) are from 
Hipparcos (ESA 1997) and Tycho-2 (tfog et al. 2000) cat- 
alogues. 

4.1. Spectral types and other derived parameters 

To obtain the spectral type of this binary system we com- 
pared our high resolution echellc spectra, in several spectral 
orders free of lines sensitive to chromospheric activity, with 
spectra of inactive reference stars of different spectral types 
and luminosity classes observed during the same observ- 
ing run. This analysis makes use of the program STARMOD 
developed at Penn State University (Barden 1985) which 
we modified later. This program constructs a synthesized 
stellar spectrum from artificially rotationally broadened, 
radial-velocity shifted, and weighted spectra of appropri- 
ate reference stars. 

For FF UMa we obtained the best fit between observed 
and synthetic spectra using a K1IV reference star for pri- 
mary component and a K0V for the secondary, with a con- 
tribution to the continuum of 0.70 and 0.30 respectively. 
These spectral types are in agreement with the results re- 
ported by other authors who suggested an evolved compo- 
nent. In our spectra, the spectral features indicate strongly 
the subgiant nature of the primary. 

We note that since the primary component is a sub- 
giant star, the stellar parameters, such us mass and ra- 
dius, cover a wide range of values. Therefore, we determined 
these crucial characteristics using data from the secondary 
star. Assuming a K0V spectral type for secondary compo- 
nent, we adopted from Landolt-Bornstein tables (Schmidt- 
Kaler 1982) a mass Ms = 0.79 M and, according to the 
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SB B-V V — R P„ rb 
(days) 


-Pphot 

(days) 


(km s" 1 ) 


7T 

(mas) 


[i a cos 8 
(mas/yr) 


fJ-s 
(mas/yr) 


KOV/KOIV 1 


2 0.97 0.75 3.274 1 


3.27 


33.57 1 ±0.45/32 .38 1 ±0.75 


9.57±0.92 


-21.20±1.30 


-23.00±1.60 



1 values determined in this paper 



mass ratio from the orbital solution (q = 2.12), we derived Table 6. Galactic space-velocity components 
a primary mass of Mp = 1.67 Mq. 

In addition, from the photometric period (3.27 days) 
given by Henry et al. (1995) and the rotational velocity, 
calculated here, i>sinip = 33.57 km s _1 (Sect. 4.2), we 
estimated a minimum radius of i?sinip — 2.17 Rq. This 
agrees with the subgiant radii and previous estimations. 



U±au 


V±a v 


W ± aw 


VTotal 


(km s" 1 ) 


(km s" 1 ) 


(km s" 1 ) 


(km s" 1 ) 


-7.04±1.14 


-12.92±1.33 


-6.68±0.68 


16.16 



4.2. Rotational velocities 



Jeffries et al. (1995) estimated the projected rotational ve- 
locity (usini) as 41 km s _1 for the primary component and 
as 32 km s _1 for the secondary. Fekel (1997) obtained 38.8 
and 39.7 km s _1 for each component and Strassmeier et al. 
(2000) reported lower values, 17 and 16 km s" 1 respectively. 

By using the program STARMOD we obtained the best 
fits for each observing run using v sin i values os ~35 km s _1 
for primary component and «38 km s~ 1 for secondary com- 
ponent. 

To determine a more accurate rotational velocity of this 
star we made use of the cross-correlation technique in our 
high resolution echelle spectra by using the routine FXCOR 
in IRAF. The method is described carefully in previous 
papers (see Galvez et al. 2002; and Lopez-Santiago et al. 
2003) and is based on the fact that when a stellar spec- 
trum with rotationally broadened lines is cross-correlated 
against a narrow-lined spectrum, the width of the cross- 
correlation function (CCF) is sensitive to the amount of 
rotational broadening of the first spectrum. 

As a template star in this process we used the K1V star 
HD 26965 for the primary component and the K0V star 
HD 3651 for the secondary for the McDonald run and the 
K2V star HD 166620 for both components in the remain- 
ing runs. All these stars have very low rotation velocity -less 
than 3 km s . The averaged values obtained are wsini = 
33.57±0.45 and 32.38±0.75 km s -1 for primary and sec- 
ondary components respectively. 




4.9x10* 4.95x10* 5x10* 5.05x10* 5.1x10* 5.15x10* 5.2x10* 5.25x10* 5.3x10* 
HJD 2400000.4- 

Fig. 5. The (O - C) (Observed - Calculated T conj ) versus 
Heliocentric Julian date HJD for every observing run. 




4.3. Kinematics 

Computing the galactic space-velocity components (U, V, 
W) of FF UMa required both radial velocity and precise 
proper motions and parallax. For the former, we used the 
center of mass velocity, 7, determined in the orbital solution 
for the FOCES04 observing run (see Sect. 3.4). For the 
latter we utilized data taken from Hipparcos (ESA 1997) 
and Tycho-2 (H0g et al. 2000) catalogues (see Table EJ) (for 
details see Montes et al. 2001a, b). 

In addition, we included FF UMa in an extended study 
of binary star kinematics in young moving groups. It in- 
cluded the application of Eggen's peculiar velocity and ra- 
dial velocity criteria (see Montes et al. 2001a and reference 
therein) and spectroscopic criteria (see the Li 1 A6707.8 line 
Sect. 4.4). 



Fig. 6. Center of mass velocity, 7, obtained in each observ- 
ing run fit versus time. 



The resulting values of (U, V, W) and associated errors 
are given in Table [H These errors have been calculated 
assuming the value of 7 determined from the FOCES04 
observing run. However, taking into account the changes in 
7 ss 3 km s^ 1 we detected, we expect the uncertainties to 
be larger (see Sect. 3.3). 

Using the (U, V) and (V, W) diagrams (Eggen 1984, 
1989; Montes et al. 2001a), the velocity components lie 
clearly inside the Castor moving group boundaries. In ad- 
dition, Eggen's radial velocity criteria also confirm their 
membership of FF UMa to this group (Galvez 2005). 
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4.4. The Lil X6707.8 line 

As it is well known, Li I A6707.8 spectroscopic feature is 
an important diagnostic of age in late-type stars, since it is 
destroyed easily by thermonuclear reactions in the stellar 
interior. 

The spectral region of the resonance doublet of Li I 
at A6708 A is covered by most of our spectral observa- 
tions. Despite blending with photospheric lines, mainly Fe I 
(6707.4 A), we could separate the contribution from both 
components. We then measured the equivalent width (EW 
hereafter) of (Li I + Fe i) of both components in our ob- 
served spectra. We calculated the contribution of Fe I by 
using both calibrations of Fe 1-effective temperature from 
Soderblom et al. (1990) and Fe i-(B — V) color index from 
Favata et al. (1993). We obtained the corrected EW {hi i) 
by subtracting the EW(Fe i) of the total measured equiv- 
alent width, EW (hi i+Fe i). The resulting mean values 
of EW{hi i) are 200 mA for the primary component and 
141 mA for the secondary. These values are corrected by 
the contribution of each component to the continuum (see 
Sect. 5). 

By using the spectral subtraction technique, that is, ob- 
taining the EW of Li I directly from the subtracted spectra, 
we obtained a mean EW(U i) of 132 mA and 86 mA for 
both component respectively. 

The discrepancy between the two methods used to cal- 
culate the EW is due to the influence of stellar metallici- 
ties. In the first technique, stellar metallicity is not taken 
into account in the relation calibrations, and in the second 
technique, the stellar metallicity of the standard star used 
to create the synthetic spectra is undetermined. In spite of 
this, these EW (hi i)s values are of the same order as the 
Li I EWs of other Castor moving group members, which 
have an age around 200 Myr. 



5. Chromospheric activity indicators 

The echelle spectra analyzed allow us to study the behav- 
ior of the Ca II H & K to the Ca n IRT lines, different 
indicators formed at varying atmospheric heights. We de- 
termined the chromospheric contribution of these features 
using the spectral subtraction technique described in detail 
by Montes et al. (1995) and Papers I, II, III and IV. We 
constructed the synthesized spectrum using the program 

STARMOD. 

Taking into account the stellar parameters derived in 
Sect. 4 we used reference stars of the K1IV spectral type 
for the primary component and K0V spectral type for the 
secondary component (see Sect. 4.3), with a contribution 
of 0.70/0.30 respectively. 

In Table[7] (only available in electronic form) we present 
the excess emission EW, measured in the subtracted spec- 
tra, for the Ca II H & K, He, RS, H7, H/3, Ha, and Ca II 
IRT lines in all observing runs. We list the EWs of emis- 
sion features for both components (P/S); when lines were 
blended we list only total EWs. The uncertainties in the 
measured EW were estimated taking into account: a) the 
typical internal precisions of STARMOD (0.5 - 2 km s _1 in 
velocity shifts, and ±5 km s _1 in vsmi), b) the rms ob- 
tained in the fit between observed and synthesized spectra 
in the spectral regions outside the chromospheric features 
(typically in the range 0.01-0.03) and c) the standard devi- 



E 



Fig. 7. Spectra in the Ha line region in McDonald98 ob- 
serving run. We plot the observed spectrum (solid-line) and 
the synthesized spectrum (dashed-line) in the left panel, 
and the subtracted spectrum (dotted line), in the right 
panel. We mark position of primary component lines with a 
(P) and position of secondary component lines with a (S). 



ations of the EW measurements. The final estimated errors 
are in the range 10-20%. 

We corrected the measured EWs for the relative contri- 
bution of each component to the total continuum (Sp and 
5s), using the radii assumed in Sect. 4.1. and temperatures 
from Landolt-Bornstcin tables (Schmidt-Kaler 1982). We 
obtained the final EWs for the components multiplying by 
a factor 1/Sp and 1/Ss, respectively. We present the result 
in Table [71 

These adopted EWs were transformed to absolute sur- 
face fluxes using the empirical stellar flux scales calibrated 
by Hall (1996) as a function of the star color index. We 
used the B — V index and the corresponding coefficients for 
Ca 11 H & K, Ha and Ca 11 IRT. We used for He the same 
coefficients as for Ca 11 H & K, and derived the US, H7 and 
H/3 coefficients of flux by making an interpolation between 
the values of Ca 11 H & K and Ha. We present the loga- 
rithm of the obtained absolute flux at the stellar surface 
(logFg) for the different chromospheric activity indicators 
in Table [5] (only available in electronic form) . 

In Figs. [7]- (HI we plot the Ha, Ca 11 H & K and Ca 11 
IRT A8498, A8542 lines region for each observation, the ob- 
served spectrum (solid-line) and the synthesized spectrum 
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Fig. 8. The same as in previous Figure in HETOO observing Fig. 9. The same as in previous Figure in FOCES02 ob- 
run. serving run. 



(dashed-line) in the left panel and the subtracted spectrum 
(dotted line) in the right panel. We included the observing 
run of each spectrum in these figures. In Fig. [THI we plot a 
representative subtracted spectrum of FF UMa in the H/3 
line region. 



The EW average value measured in the subtracted spec- 
tra is EW(Ha) — 1.64/2.47 A for the primary and sec- 
ondary components. We note that these are higher values 
than those reported by Jeffries et al. (1995). In Table[7jwe 
list the EW of each stellar component determined by the 
fit described above. We also list the total EW (primary + 
secondary) determined by integrating the total excess emis- 
sion profile. We note that Ha line shows notable variations 
with orbital phase but also from one epoch to another in 
both components. 



5.1. Ha 

The Ha line is observed in emission above the continuum in 
all the spectra (see Figs. [71- [Til left panel). In the observed 
spectrum, the emission associated with the secondary is 
larger than that associated with the primary. However, after 
applying the spectral subtraction technique, the Ha emis- 
sions above the continuum coming from both components 
are similar and in some cases the primary one is larger. The 
Ha emission is persistent during all observations indicating 
that it is a very active binary system similar to RS CVn 
and BY Dra systems that always show Ha emission above 
the continuum. Measuring the EW of this line, we found 
that each stellar component is formed by a central nar- 
row component and a broad component that moves from 
red to blue. These are an indication of microflare activity 
(see Papers I, II, and III). While we were able to sepa- 
rate the narrow components, we were unable to deblend 
the broad ones. Therefore, to determine the contribution of 
each stellar component to the total excess emission, we fit 
the narrow and broad component of each star together (see 
Figs. Ql and Ql|). 



5.2. H/3, H 7 and HS 

We can see absorption of H/3, H7 and US Balmer lines 
filled in with emission in the observed spectra. After ap- 
plying the spectral subtraction, clear excess emission is de- 
tected from both components (see a representative spec- 
trum in the H/3 line region in Fig. [T5)) . When the S/N 
was high enough we deblended the emission coming from 
both components by using a two-Gaussian fit to the sub- 
tracted spectra (see Table [7J . These three lines show the 
same behavior with orbital phase that the Ha line in both 
components. Their mean values are EW(Hf3) — 0.33/0.25 
A, EW(Hj) = 0.14/0.19 A and EW(RS) = 0.13/0.15 A. 
We also measured the ratio of excess emission in the Ha 



and H/3 lines, ^^§^y 7 and the ratio of excess emission 



with the correction: 



given by Hall & Ramsey (1992). This takes into account the 
absolute flux density in these lines and the color difference 
in the components. We obtained mean values of ~ 
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Fig. 10. The same as in previous Figure in FOCES04 ob- 
serving run. 

6 for the primary component and w 5 for the secondary. 
These values indicate, according to Buzasi (1989) and Hall 
& Ramsey (1992), the presence of prominence-like material 
above the stellar surface in both components of the system. 

5.3. Can H & K and He 

The Ca II H k K line region is included in FOCES 2002 
and 2004 and NOT04 observing runs. 

This spectral region is located at the end of the echel- 
logram, where the efficiency of the spectrograph and the 
CCD decrease very rapidly and therefore the S/N ratio ob- 
tained is very low; thus the normalization of the spectra is 
very difficult. In spite of this, the spectra show strong emis- 
sion in the Ca II H & K lines and a clear emission in the 
He line from both components (see Fig. [T2"|) . These allow us 
to apply the spectral subtraction in this spectral region. As 
we can see in Fig. [T^l the He line arising from one of the 
component overlaps with the Ca n H line arising from the 
other component at some orbital phases, so their EW were 
measured with a Gaussian fit when it was possible. Mean 
EWs values measured in these spectra are EW= 1.23/1.14 
A for each component in Ca II K line, EW— 1.23/1.08 A 
in Ca II H line and EW= 0.30/0.43 A in He line. 

As the Ha emission line, the Ca II H & K lines show 
variations with both orbital phase and from one epoch to 
another in both components. 
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Fig. 11. The same as in previous Figure in NOT04 observ- 
ing run. 



5.4. Ca II IRT lines (X8498, X8542 and X8662) 

All our echelle spectra include the three lines of the Ca II 
infrared triplet (IRT) except for the A8498 line in HET00 
and NOT04 runs. In all of the spectra we observed a clear 
emission above the continuum in the core of the Ca II IRT 
absorption lines (see Figs. [T3l- fl~5|) from both components. 
After applying the spectral subtraction, we could see that 
the emission coming from the primary component is larger 
than the emission from the secondary. 

We measured mean EWs for these three Ca n 
lines of £W(A8498, A8542, A8662)= 0.58/0.35, 0.67/0.43 
and 0.59/0.35 A respectively. For each component we found 
considerable variations with orbital phase that appear anti- 
correlated with the variations in the Balmer lines. 

In addition, we calculated the ratio of excess emission 
EW. ~^ 8542 , which is also an indicator of the type of chro- 

£,8498 ' . 

mospheric structure that produces the observed emission; 
in solar plages, values of saMa ~ 1.5-3 are measured, while 
in solar prominences the values are « 9, the limit of an 
optically thin emitting plasma (Chester 1991). We found 
for this star small values of the ^ 8542 ratio, « 1.0, for both 

£,8498 

components (see Table [8]). This indicates that the Ca II 
IRT emission of this star arises from plage-like regions at 
the stellar surface, in contrast with the Balmer lines that 
come from prominences. This markedly different behavior 
of the Ca II IRT emission has also been found in other 
chromospherically active binaries (see Papers III; IV and 
references therein). 
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Fig. 13. The same as in previous figures but in Ca n IRT 
(A8498 & A8542) line regions in McDonald98 observing run. 



5.5. Variation of activity with time 

As we mentioned above, the EW emission lines that arc 
chromospheric indicators show variations with orbital phase 
due to activity features present in both stellar surfaces. But 
there are also variations from one epoch to another in both 
components. To study if there is a correlation between ac- 
tive cycle and orbital period variation as Lanza (2006) sug- 
gested for the case of HR 1099, we must have a follow up of Fig. 15. The same as previous Figure in FOCES04 observ- 
this system during a complete cycle of its variation period ing run. 
(> 22 years). This kind of study would provide us another 
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Fig. 16. Example of the Ha region fit in the subtracted 
spectrum of the two components by using IRAF SPLOT 
task. 
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Fig. 17. Another example of Ha region fit. 
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Fig. 18. The same as Ha figures but in H/3 line region in 
FOCES04 observing run. 



clue for understanding and testing the Applegate's mecha- 
nism. 



6. Conclusions 

In this paper, we present a detailed spectroscopic analysis of 
a X-ray /EUV selected chromospherically active binary sys- 



tem 2RE J0933+624 (FF UMa). We analyzed high resolu- 
tion echelle spectra that include the optical chromospheric 
activity indicators from the Ca II H & K to Ca n IRT lines, 
as well as the Li I A6707.8 line and other photospheric lines 
of interest. 

With a large number of radial velocities from the litera- 
ture and from our spectra taken over several years, we found 
that this system shows an orbital period variation similar 
to those previously found in other RS CVn systems. 

Although the existence of an unseen distant third star 
as a component of the system cannot be completely ruled 
out as the cause of these variations, we think that the 
more plausible explanation is Applegate's mechanism, or 
at least the qualitative idea; which states that the orbital 
period change is due to the gravitational coupling of the 
orbit to changes in the quadrupole moment of the magneti- 
cally active stellar components of the system. In the case of 
FF UMa, we calculated a (O — C) (T con j) that gives us a rel- 
ative orbital period variation of dP/P « 1CP 4 in 11 years, 
that is, one order of magnitude higher than the variations 
in HR 1099, the largest observed until now. We suggest here 
that this order of magnitude difference between the period 
variations in FF UMa and HR 1099 could be explained by 
the different activity level. The components of FF UMa are 
very active and have more effective dynamo mechanisms 
than HR 1099 components. 

Once we adopted an orbital period, from the FOCES04 
observing run, we improved the determination of the orbital 
solution of the system relative to previous determinations 
by other authors. We obtained a nearly circular orbit with 
an orbital period very close to photometric period, indicat- 
ing that it has a synchronous rotation. 

The spectral classifications derived by comparing FF 
UMa with spectra of reference stars, leads us to consider 
the primary component as a subgiant star and the sec- 
ondary component as a KOV star. The results from orbital 
parameters and photometric characteristics help us to ob- 
tain physical parameters from the primary, Mp = 1.67 M 
and i?sinip = 2.17 Rq, in agreement with subgiant radii 
and previous estimates. 

By using the information provided by the width of the 
cross-correlation function we determined a projected rota- 
tional velocity, v sin i, of 33.57±0.45 km s _1 and 32.38±0.75 
km s _1 for the primary and secondary components respec- 
tively. 

The presence of the Li I line is in agreement with the 
kinematics results, i.e., it belongs to the young disk and is 
probably a member of the Castor moving group. 

The study of the optical chromospheric activity indi- 
cators shows that FF UMa system has a high level of ac- 
tivity in both components. The variation of Ha and the 
rest of the Balmer and Ca II H&K lines are very similar 
and anti-correlated in phase with Ca n IRT emission, as we 
confirmed with the -f^ and f 8542 results. This indicates 

^8498 

that the Balmer emission lines arise from prominence-like 
material while the emission of Ca n IRT lines arise from 
plage-like regions. In addition, both components show vari- 
ations from one epoch to another that could have a correla- 
tion with the orbital period variation. Future spectroscopic 
and photometric studies of this system could confirm this 
hypothesis and provide a better understanding and test of 
Applegate's mechanism. 
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Table 7. EW of chromospheric activity indicators 



EW(A) in the subtracted spectra 



Obs. <p Ca II Ca II IRT 



Idt. 5 




K 


H 


He 


H<5 


H7 


H/3 


Ha 
Ha,(P+S) 


A8498 


A8542 


A8662 


(1) 


0.75 


_ 


_ 


_ 


_ 


_ 


_ 


0.89/2.45 


0.35/0.23 


0.53/0.41 


0.48/0.28 


















2.66 








(1) 


0.06 














3.64 1 


0.63 1 


1.02 1 


0.77 1 
















3.27 








(1) 


0.37 














1.54/2.93 


0.58 1 


1.18 1 


0.55 1 


















3.31 








(1) 


0.67 














1.24/2.02 
























2.43 








(1) 


0.98 














3.52 1 


0.54 1 


0.88 1 


0.66 1 


















2.96 








(1) 


0.30 


_ 






_ 






1.43/2.21 
2.88 


0.38/0.22 


0.60/0.43 


0.48/0.27 


(1) 


0.60 








_ 






0.90/2.31 


0.43/0.10 


0.80 1 


0.53/0.28 


















2.47 








(1) 


0.90 














1.00/2.15 


3 


3 


3 


















2.62 








(1) 


0.21 


_ 


_ 


_ 


_ 


_ 


_ 


1.25/2.19 


0.37/0.23 


0.70/0.56 


0.42/0.31 


















2.70 








(1) 


0.55 


_ 


_ 




_ 


_ 


_ 


1.25/1.93 


0.57 1 


0.94 1 


0.53 1 


















3.07 








(2) 


0.74 


- 


- 


- 


- 


- 


0.31/0.23 


1.95/2.73 


- 


0.71/0.32 


0.58/0.32 


















4.90 








\ Z ) 


u.u / 












u.oy 


3 42 1 




1 

L.XJ i 


u.yo 


















2.93 








(2) 


0.47 


- 


- 


- 


- 


- 


0.65 1 


3.53 1 


- 


1.16 1 


1.03 1 


















3.38 








(2) 


0.30 


_ 


_ 


- 


- 


- 


0.88 1 


1.62/3.09 


- 


0.74/0.41 


0.72/0.14 


















3.23 








(2) 


0.26 


- 


- 


- 


- 


- 


1.13 1 


5.20 1 


- 


1.25 1 


1.14 1 
















4.03 








(2) 


0.62 


- 


- 




- 


- 


0.68 1 


3.95 1 




l.ll 1 


0.95 1 


















3.33 








(2) 


0.91 


- 


- 


- 


- 


- 


0.66 1 


2.58/1.58 


- 


0.68/0.30 


0.59/0.43 


















3.19 








(2) 


0.26 


- 


- 


- 




- 


0.91 1 


3.42/1.57 
3.45 


- 


1.26 1 


1.34 1 


(3) 


0.56 


1.80/1.96 


1.46/2.17 


3 


3 


0.07/0. 13 2 


0.46/0.34 


2.15/1.99 


0.80/0.60 


0.55/0.43 


0.67/0.42 


















3.70 








(3) 


0.86 


3.42 1 


2.45 1 


0.62 1 


0.20 1 


0.26 1 ' 2 


0.74 1 ' 2 


3.99 1 


1.20 1 


0.93 1 


1.04 1 
















3.34 








(3) 


0.17 


1.92/1.27 


1.31/0.98 


0.589/ 3 


0.09/0.10 


0.12/0. 15 2 


0.45/0.21 2 


1.33/2.01 


0.91/0.38 


0.63/0.22 


0.73/0.32 


















2.91 








(3) 


0.50 


1.88/1.81 


1.53/1.11 


0.41/0.07 


0.10/0.11 


0.07/0.16 2 


0.21/0.45 2 


1.85/1.81 


0.78/0.65 


0.49/0.39 


0.73/0.43 


















3.25 








(4) 


0.70 


1.09/0.82 


1.08 4 /0.71 


0.51/1. 08 4 


0.13/0.14 


0.15/0.15 


0.36/0.34 


1.34/2.20 


0.57/0.25 


0.75/0.42 


0.66/0.29 


















2.81 








(4) 


0.26 


0.92/0.85 


0.95/0.92 4 


0.92 4 /1.49 


3 


0.18/0.08 


0.31/0.34 


1.89/1.60 


0.55/0.41 


0.70/0.44 


0.63/0.40 


















2.79 








(4) 


0.57 


3.03 1 


2.07 1 


0.24 1 


0.69 1 


0.50 1 


0.67 1 


3.45 1 


0.91 1 


1.19 1 


1.08 1 


















2.70 








(4) 


0.62 


1.40/0.71 


0.94 4 /0.97 


0.57/0. 94 4 


3 


3 /0.19 


0.36/0.28 


1.19/2.26 


0.57/0.28 


0.60/0.41 


0.56/0.43 


















2.58 








(4) 


0.66 


3 


3 


3 


3 


0.15/0.21 


0.30/0.27 


1.19/2.06 


0.60/0.31 


0.72/0.46 


0.58/0.36 


















2.59 








(4) 


0.87 


3.12 1 


2.14 1 


0.21 1 


0.31 1 


0.37 1 


0.71 1 


3.47 1 


0.85 1 


1.25 1 


1.05 1 
















2.78 








(4) 


0.91 


3.04 1 


1.90 1 


0.33 1 


0.35 1 


0.24 1 


0.60 1 


3.31 1 


0.86 1 


1.16 1 


0.95 1 
















2.97 








(4) 


0.95 


3 


1.34 1 


0.31 1 


3 


0.20 1 


0.63 1 


3.46 1 


0.90 1 


1.15 1 


0.95 1 


















3.03 








(4) 


0.18 


0.91/0.76 


1.09/0.98 4 


0.98 4 /1.12 


3 


0.16/0.16 


0.22/0.31 


1.86/1.67 


0.58/0.40 


0.70/0.49 


0.68/0.38 


















3.07 








(4) 


0.22 


1.52/1.37 


1.33/0.99 4 


0.99 4 /1.70 


0.16/0.24 


0.17/0.29 


0.27/0.31 


2.09/1.59 


0.59/0.41 


0.80/0.53 


0.69/0.43 


















3.19 








(4) 


0.27 


1.25/1.24 


1.08/0.96 4 


0.96 4 /0.36 


0.17/0.15 


0.18/0.23 


0.31/0.30 


1.99/1.74 


0.57/0.42 


0.75/0.53 


0.64/0.43 


















3.10 








(4) 


0.48 


3.30 1 


1.97 1 


0.44 1 


0.50 1 


0.29 1 


0.67 1 


3.41 1 


0.94 1 


1.28 1 


1.12 1 
















2.91 








(4) 


0.56 


2.67 1 


2.22 1 


0.28 1 


0.36 1 


0.44 1 


0.70 1 


3.55 1 


0.99 1 


1.36 1 


1.18 1 
















2.91 








(5) 


0.29 


1.33/1.20 


0.77/1.66 4 


1.66 4 


3 


3 


0.26/0.26 


2.04/1.56 




0.91/0.43 


0.34/0.46 


















3.16 








(5) 


0.27 


1.16/1.02 


0.54/0.66 4 


0.66 4 


3 


3 


0.36/0.28 


1.88/0.59 




0.60/0.64 


0.51/0.40 


















2.87 








(5) 


0.64 


1.09/0.59 


0.89 4 /0.81 


4 /0.35 


3 


3 


0.30/0.24 


2.05 




0.65/0.26 


0.52/0.22 


















2.87 








(5) 


0.21 


0.18/1.00 


1.10/0.96 4 


4 /0.48 


0.51 1 


0.14/0.11 


0.41/0.21 


1.89/1.83 




0.59/0.42 


0.58/0.40 


















3.19 
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Table 8. Emission fluxes 



logFs 


Obs. 




Ca II 












Ca II IRT 




Idt. 5 


K 


H 


He 


H<5 


H 7 H/3 


Ha 


A8498 


A8542 


A8662 



(!) 














6.89/7.28 


6.41/6.13 


6.59/6.38 


6.55/6.22 


(1) 














7.50 1 


6.67 1 


6.88 1 


6.76 1 


(1) 


_ 


_ 


_ 




_ 


_ 


7.13/7.36 


6.63/6.54 


6.94/6.84 


6.61/6.51 


(1) 


- 


- 


- 


- 


- 


- 


7.04/7.20 
7.49 1 


- 


- 


- 


(1) 


- 


- 


- 


- 


- 


- 


6.60 1 


6.81 1 


6.69 1 


(1) 






- 


- 






7.10/7.24 


6.45/6.11 


6.65/6.41 


6.55/6.20 


(1) 


- 


- 


- 


- 


- 


- 


6.90/7.25 


6.50/6.77 


6.77/6.67 


6.59/6.22 


(1) 


- 


- 


- 


- 


- 


- 


6.94/7.22 


3 


3 


3 


(1) 


- 


- 


- 


- 


- 


- 


7.04/7.23 


6.44/6.13 
6.63 1 


6.72/6.52 


6.49/6.26 


(1) 


- 




- 


- 


- 




7.04/7.58 


6.84 1 


6.59 1 


(2) 


- 




- 


- 




6.40/6.20 
6.68 1 


7.23/7.33 




6.72/6.28 
6.90 1 


6.63/6.28 


(2) 


- 




- 


- 


- 


7.48 1 


- 


6.84 1 


(2) 


- 


- 


- 


- 


- 


6.85 1 


7.15/7.38 


- 


6.74/6.38 
6.93 1 


6.73/5.92 


(2) 


- 


- 


- 


- 


- 


6.72 1 


7.49 1 


- 


6.88 1 


(2) 


- 


- 


- 


- 


- 


6.74 1 


7.54 1 


- 


6.92 1 


6.85 1 


(2) 


- 


- 


- 


- 


- 


6.96 1 


7.66 1 


- 


6.97 1 


6.93 1 


(2) 


- 


- 


- 


- 


- 


6.73 1 


7.35/7.09 


- 


6.70/6.25 


6.64/6.41 


(2) 


- 


- 


- 


- 


- 


6.86 1 


7.48/7.09 


- 


6.97/6.87 


7.00/6.90 


(3) 


7.15/7.09 
7.43 1 


7.06/7.13 
7.28 1 


3 


3 


5.73/5.94 
6.30 1 


6.57/6.38 
6.77 1 


7.28/7.19 
7.54 1 


6.77/6.55 
6.9 5 1 


6.61/6 .40 
6.84 1 


6.70/6.40 
6.89 1 


(3) 


6.69 1 


6.17 1 


(3) 


7.18/6.90 


7.01/6.79 


6.67/ 3 


5.83/5.82 


5.96/6.00 


6.56/6.16 


7.07/7.19 


6.83/6.35 


6.67/6.11 


6.73/6.28 


(3) 


7.17/7.06 


7.08/6.84 


6.51/5.64 


5.87/5.86 


5.73/6.03 


6.23/6.49 


7.21/7.15 


6.76/6.58 


6.56/6.36 


6.73/6.41 


(4) 


6.93/6.71 


6.93/6.65 


6.60/ 3 


5.99/5.96 


6.06/6.00 


6.46/6.37 


7.07/7.23 


6.63/6.17 


6.75/6.40 


6.69/6.23 


(4) 


6.86/6.73 


6.87/6.76 
7.21 1 


3 /6.97 




6.14/5.73 
6.58 1 


6.40/6.37 
6.73 1 


7.22/7.10 

7.48 1 


6.61/6.38 
6.83 1 


6 .72/6.42 
6.95 1 


6.67/6.37 
6.90 1 


(4) 


7.38 1 


6.28 1 


6.71 1 


(4) 


7.04/6.65 


6.872/6.79 


6.65/ 3 


3 


3 /6.10 


6.46/6.28 


7.02/7.25 


6.63/6.22 


6.65/6.38 


6.62/6.41 


(4) 


3 


3 


3 


3 


6.06/6.15 


6.38/6.27 


7.02/7.20 


6.65/6.26 


6.73/6.43 


6.63/6.33 


(4) 


7.39 1 


7.23 1 


6.22 1 


6.37 1 


6.45 1 


6.76 1 


7.48 1 


6.80 1 


6.97 1 


6.89 1 


(4) 


7.38 1 


7.17 1 


6.41 1 


6.42 1 


6.26 1 


6.68 1 


7.46 1 


6.80 1 


6.93 1 


6.85 1 


(4) 


3 


7.02 1 


6.39 1 


3 


6.18 1 


6.70 1 


7.48 1 


6.82 1 


6.93 1 


6.85 1 


(4) 


6.85/6.68 


6.93/6.79 4 


3 /6.85 


3 


6.09/6.03 


6.25/6.33 


7.21/7.11 


6.63/6.37 


6.72/6.46 


6.70/6.35 


(4) 


7.08/6.94 


7.02/6.79 4 


3 /7.03 


6.08/6.20 


6.11/6.29 


6.34/6.33 


7.26/7.09 


6.64/6.38 


6.77/6.50 


6.71/6.41 


(4) 


6.99/6.89 


6.93/6.78 4 


3 /6.35 


6.10/5.99 


6.14/6.18 


6.40/6.31 


7.24/7.13 


6.62/6.40 


6.75/6.50 


6.68/6.41 


(4) 


7.41 1 


7.19 1 


6.54 1 


6.57 1 


6.35 1 


6.73 1 


7.48 1 


6.84 1 


6.98 1 


6.92 1 


(4) 


7.32 1 


7.24 1 


6.34 1 


6.43 1 


6.53 1 


6.75 1 


7.49 1 


6.87 1 


7.00 1 


6.94 1 


(5) 


6.78/7.02 


7.02/6.88 4 


7.12/7.02 


3 /3 


3 /3 


6.32/6.25 


7.25/7.08 




6.83/6.41 


6.40/6.43 


(5) 


6.63/6.62 


6.96/6.81 4 


6.71/ 3 


3 /3 


3 /3 


6.46/6.28 


7.22/6.66 




6.65/6.58 


6.58/6.37 


(5) 


6.84/6.71 


6.93 4 /6.57 


4 /6.34 


3/3 


3/3 


6.38/6.21 


7.00/7.20 




6.68/6.19 


6.59/6.11 


(5) 


6.94/6.78 


6.15/6.80 4 


6.90/6.48 


6.58 1 


6.03/5.86 


6.52/6.16 


7.22/7.15 




6.49/6.40 


6.63/6.37 



Notes as in previous Table. 



